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Abstract 

In this paper, we study resource allocation and adaptive modulation in SC-FDMA which is adopted as 
the multiple access scheme for the uplink in the 3GPP-LTE standard. A sum-utility maximization (SUmax), 
and a joint adaptive modulation and sum-cost minimization (JAMSCmin) problems are considered. Unlike 
OFDMA, in addition to the restriction of allocating a sub-channel to one user at most, the multiple 
sub-channels allocated to a user in SC-FDMA should be consecutive as well. This renders the resource 
allocation problem prohibitively difficult and the standard optimization tools (e.g., Lagrange dual approach 
widely used for OFDMA, etc.) can not help towards its optimal solution. We propose a novel optimization 
framework for the solution of these problems that is inspired from the recently developed canonical duality 
theory. We first formulate the optimization problems as binary-integer programming problems and then 
transform these binary-integer programming problems into continuous space canonical dual problems that 
are concave maximization problems. Based on the solution of the continuous space dual problems, we derive 
resource allocation (joint with adaptive modulation for JAMSCmin) algorithms for both the problems which 
have polynomial complexities. We provide conditions under which the proposed algorithms are optimal. 
We also propose an adaptive modulation scheme for SUmax problem. We compare the proposed algorithms 
with the existing algorithms in the literature to assess their performance. 

The authors are with the Department of Telecommunications, Ecole Superieure d'Electricite (Supelec), 91192 Gif-sur-Yvette, 
France (e-mail: {ayaz.ahmad, mohamad.assaad}@supelec.fr). 
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I. Introduction 

Single Carrier Frequency Division Multiple Access (SC-FDMA) is currently attracting a lot of attention 
as an alternative to OFDMA in the uplink. Its low PAPR feature has the potential to benefit the mobile 
terminals in term of transmit power efficiency. In fact, SC-FDMA is a single carrier multiple access 
technique which utilizes single carrier modulation and frequency domain equalization. Its overall structure 
and performance are similar to that of OFDMA system. Unlike the parallel transmission of the orthogonal 
sub-channels in OFDMA, the sub-channels are transmitted sequentially in SC-FDMA. This sequential 
transmission of sub-channels considerably reduces the envelope fluctuation in transmitted waveform and 
results in low PAPR ID. There are two types of SC-FDMA: localized-FDMA (L-FDMA) in which the 
sub-channels assigned to a user are adjacent to each other, and interleaved-FDMA (I-FDMA) in which 
users are assigned with sub-channels distributed over the entire frequency band [1]. In 3GPP-LTE standard 
0, the current working assumption is to use OFDMA for downlink and localized SC-FDMA for uplink. 

A. SC-FDMA vs OFDMA from a Resource Allocation Perspective 

Most of the previous work on resource allocation has focused on power and sub-channels allocation in 
downlink OFDMA systems 0- Q. One of the well known approaches for solving the OFDMA resource 
allocation problem is exploiting its time-sharing property |8]. Based on this property, it is shown in (7), 
and El that for practical number of sub-channels, the resource allocation problem in OFDMA systems 
can be solved by Lagrange multipliers method with zero duality gap. However, none of above is directly 
applicable to uplink SC-FDMA. This is due to the fact that in localized SC-FDMA in addition to the 
restriction of allocating a sub-channel to one user at most, the multiple sub-channels allocated to a user 
should be adjacent to each other as well. Furthermore, a frequency domain equalizer is used in SC-FDMA 
over all the sub-channels allocated to the user which makes the SNR expression much more complicated 
than in OFDMA where the SNR on each sub-channel is independent from the other sub-channels. 

The common approach used for resource allocation in OFDMA is to formulate the mutual exclusivity 
restriction on sub-channels allocation as binary-integer constraint, solve the problem to get an approximated 
solution in continuous domain, and then discretize the continuous values into the closest binary values. 
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But in SC-FDMA resource allocation, this approach cannot be employed. The reason is that if the problem 
is solved by relaxing the 0-1 constraint, then, during discretization of the continuous domain solution, 
the adjacency constraint on sub-channels allocation cannot be assured. This necessitates the design of a 
framework that also ensures the adjacency constraint on sub -channels allocation which is a very difficult 
task. 

B. Related Work 

In most of the previous work on SC-FDMA, the implementation problems in the physical layer are 
studied (e.g., 0- iTTTlD . The resource allocation problem in uplink SC-FDMA has also been addressed 
in a number of publications. In |[T2l . a heuristic opportunistic scheduler for allocating frequency bands 
to the users in the uplink of 3G LTE systems is proposed. In [fT3l, the authors have proposed a greedy 
sub-optimal schedular for uplink SC-FDMA systems that is based on marginal capacity maximization. In 
|[T4l . the authors revise the same framework used in ||T3l for developing a proportional fair scheduling 
scheme. However, in addition to being sub-optimal, the proposed schedulers in both ||T3l and |[T4l do 
not consider the sub-channels adjacency constraint which is an important physical layer requirement for 
localized SC-FDMA. In |fI51l , a set of greedy sub-optimal proportional fair algorithms for localized SC- 
FDMA systems is proposed in the frequency-domain setting. This work respects the sub-channels adjacency 
constraint but does not consider any constraint on the power. In [16], a weighted-sum rate maximization 
in localized SC-FDMA systems is considered where the problem is formulated as a pure binary-integer 
program. Though the proposed binary-integer programming framework captures all the basic constraints of 
the localized SC-FDMA and allows to perform resource allocation without resorting to exhaustive search, 
it is still not the best solution as the 0-1 requirement turns the problem into combinatorial with exponential 
complexity. Thus, keeping in view the computational complexity of the binary-integer programming, the 
authors have also proposed a greedy sub-optimal algorithm that is similar in spirit to the approach in ||T3l 
with an additional constraint on the adjacency of the allocated sub-channels. Moreover, all the previous 
work is based on rate/capacity maximization and no work to the best of our knowledge has considered 
power minimization in uplink SC-FDMA systems. Since the mobile terminals have limited energy, energy- 
economization is needed and fast power control should be considered while allocating the resources to the 
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users in the uplink. 

C. Motivation and Contributions 

In this paper, we consider resource allocation and adaptive modulation in localized SC-FDMA systems. 

We consider two optimization problems: sum-utility maximization (SUmax), and joint adaptive modulation 
and sum-cost minimization (JAMSCmin). Both these problems are combinatorial in nature whose optimal 
solutions are exponentially complex in general. The performance metric considered in the SUmax problem 
is the total utility of the system. Utility is basically an economics concept that reflects the user satisfaction 
in the system. We assume that each user in the system has an associated utility function, and the objective 
is to propose a polynomial-complexity resource allocation framework that could maximize the sum-utility 
while respecting all the constraints of localized SC-FDMA systems specific to the LTE uplink. The user 
utility function specific to this paper is defined as an arbitrary function that is monotonically increasing in 
user's SNR. The performance of the system can be further enhanced by choosing an efficient modulation 
scheme for each user. Therefore, based on the resource allocation, we also propose an adaptive modulation 
scheme, wherein an appropriate modulation is chosen for each user depending upon its effective SNR. The 
cost associated to each user in the JAMSCmin problem is a function that is monotonically increasing in 
the transmit power of that user. The objective of the JAMSCmin is to propose a low-complexity framework 
that jointly allocates the transmit powers, sub-channels and the modulation schemes to the users in order 
to minimize the total transmit power while ensuring the individual target data rates of the users as well as 
capturing the basic constraints of the localized SC-FDMA systems. The joint adaptive modulation in the 
JAMSCmin problem is important due to the fact that in order to ensure the target data rate of the users, the 
powers and sub-channels allocation should take into account the modulation schemes used by the users. 

In this paper, we propose a novel framework for the solution of the above problems. In our optimization 
framework, first we formulate the optimization problems as binary-integer programming problems. We 
then transform the binary-integer programming problems into canonical dual problems [18] in continuous 
space that are concave maximization problems under ceratin conditions. We provide the global optimality 
conditions under which the solution to each dual canonical problem is identical to the solution of the 
corresponding primal problem. We also explore some bounds on the sub-optimality of the proposed 
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framework when the optimality conditions are not satisfied. Our proposed framework has polynomial 
complexity which is a significant improvement over exponential complexity. 

The rest of this paper is organized as follows: Sections II provides the system model, and Section III 
presents the problems formulation. The canonical dual optimization framework for the solution of both the 
problems is provided in Section IV, and the resource allocation and adaptive modulation algorithms are 
derived in Section V. Section VI illustrates the numerical results, and Section VII concludes the paper. 

The following notations are used in this paper. Superscripts (.) T , and (.) H stand for transpose, and 
Hermitian of a vector or a matrix respectively. Uppercase and lowercase boldface letters denote matrices, 
and vectors respectively. The word "dual" used in this paper refers to "canonical dual". 

II. System Model 

We consider the uplink of a single cell model that utilizes localized SC-FDMA. The generalization to 
multi-cell scenario is straightforward by considering the inter cell interference in the signal-to-interference- 
plus-noise ratio (SINR) expression. We make it clear that this paper does not study inter-cell interference 
reduction/mangement but aims to optimize the resources in each cell by an efficient resource allocation 
algorithm. In the cell, K users are summed to be simultaneously active. The total bandwidth B is divided 
into N sub-channels each having 12 sub-carriers. The channel is assumed to be slowly fading or in other 
words assumed to exhibit block fading characteristics. The coherence time of the channel is greater than the 
transmission-time-interval (TTI) so that the channel stays relatively constant during the TTI (in 3GPP-LTE, 
TTI = 0.5msec). The users' channel gains are assumed to be perfectly known. 

In the following, all signals are represented by their discrete time equivalents in the complex baseband. 
Assume that Nk be the number of consecutive sub-channels allocated to user k (since a sub-channel 
cannot be allocated to more than one user simultaneously, Ylk=i = N). Let Sk = [sk,i, Sfc,7vJ T be 
the modulated symbol vector of the kth user, and Fn and denote an iV-point DFT and an N -point 
Inverse DFT (IDFT) matrices respectively. The assignment of the data modulated symbols s& to the user 
specific set of Nk sub-channels can be described by a Nk -point DFT precoding matrix FAr fc , a N * Nk 
mapping matrix and an iV-point IDFT matrix F^. The mapping matrix represents the blockwise 



sub-channel allocation where the elements D k (n,q) for n = 0, ...,N — 1 and q = 0, ...,N k — 1 are given 
by 

D k (n,q)={ ^ 3 ~ X (1) 



The transmitted signal is then 



elsewhere 



x fe = F^DfcFVSfc (2) 



At the receiver, the received signal is transformed into the frequency domain via a iV-point DFT. The 
received signal vector for user k assuming perfect sample and symbol synchronization, is given as 

y k = H fe Ff D k ¥ Nk s k + z k (3) 

where H k = diag(/ifc i, h k ^) and z k = [z kt x, z k ^Y are respectively the diagonal channel response 
matrix and the diagonal Additive White Gaussian Noise (AWGN) vector in the frequency domain. A 
frequency domain equalizer is then used in order to mitigate the ISI. The equalized symbols are transformed 
back to the time domain via an N k -point IDFT, and the detection takes place in the time domain. Let 
P k: n, and a 2 z denote the transmit power of user k on sub-channel n, and the ambient noise variance at the 
receiver for user k respectively. After several manipulations, the effective SNR for user k can be obtained 
as follows ifTTTl : 

ZF ( 1 1 ^ MMSE ( \ i\ 

lk \N k ^ P kn G kn ' 7fc I _L yN k P k , n G k . n I w 

where ^y k F is the SNR when ZF equalizer is used and is the SNR when MMSE equalizer is 

\h I 2 

used, and where G k ^ n = 1 ^ . The optimization framework proposed in this paper assumes an MMSE 
frequency domain equalization at the receiver. Nevertheless, the proposed framework is equally applicable 
for ZF equalization at the receiver. 

Unlike OFDMA where a different constellation can be adopted for each sub-channel, in SC-FDMA a 
single constellation is chosen for each user depending upon its channel quality. This is due to the fact that the 
transmit symbols directly modulate the sub-channels in OFDMA whereas in SC-FDMA, the transmit sym- 
bols are first fed to the FFT block and the output discrete Fourier terms are then mapped to the sub-channels. 
In 3GPP LTE, the constellation for each user is chosen from the set M = {QPSK, 16QAM, 64QAM}. 



III. Problems Formulation 

In this section, we formulate the two optimization problems and their equivalent binary-integer program- 
ming (BIP) problems respectively. The formulation of the problems as equivalent binary integer programs 
is an intermediate step towards its solution which are then approached by the canonical dual method. 

A. Sum-Utility Maximization (SUmax) 

1 ) SUmax Problem Formulation: We want to maximize the sum-utility subject to constraint on the total 
transmit power of each individual user p™ ax . We also have per sub-channel peak power constraint, P^ k 
i.e., the peak power transmitted on each sub-channel by any user should not exceed P^ k so that the PAPR 
is kept low |2|. In addition, in SC-FDMA for LTE uplink, the power on all the sub-channels allocated to 
a user should be equal [2], so that the low PAPR benefits could retain 12. The utility of user k denoted 
as Uki'jk) is an arbitrary function that is monotonically increasing in user's SNR The overall resource 
allocation problem can be formulated as 

K 



max 



£17* (7k) (5) 

k=l 

S-t. £ P k , n < P^™, Vfc 

n€A4 

Pk,n<Pl™ k > V ^' n 
Pk,n = Pk,h Vk,n,l 

Af k nN j = tt,Vk^j 
nfl( [j Mj) = I n G {n 1} m + l,...,n 2 - l,n 2 } >,VA; 

where A4 with cardinality N k is the set of sub-channels allocated to users k, m = min(A/fc) and rii = 
max(A/fc)- The fourth constraint determines that each sub-channel is allowed to be allocated to one user 
at most while the last constraint ensures that the sub-channels included in the set J\f k are consecutive. The 
optimization problem © is combinatorial in nature. There is a twofold difficulty in solving this problem, 
that is in addition to the exclusivity restriction on the sub-channel allocation, the allocated sub-channels 
to any user should be adjacent as well. For example, for K = 10 users and N = 24 sub-channels, the 



optimal solution requires a search across 5.26 x 10 12 possible sub-channel allocations fi~6l . which is not 
practical. 

2) Equivalent BIP Problem for SUmax problem: As an intermediate step towards its solution, we 
transform the problem to a binary-integer programming where the decisions are made on the basis of 
feasible set of sub-channel allocation patterns that satisfies the exclusivity and adjacency constraints and 
not on the basis of individual sub-channels. In other words, we form groups of contiguous sub-channels 
which will be optimally allocated among the users while respecting the exclusive sub-channels allocation 
constraint. The idea of allocation of sub-channel patterns is the same as in |16|. We elaborate the general 
idea of forming the feasible sub-channel patterns with a small example. Let us suppose that we have K = 2 
users and N = 4 sub-channels. In any allocation pattern, we put 1 if a sub-channel is allocated to a user, 
and put if it is not allocated to the user. Thus, keeping in view the sub-channel adjacency constraint, the 
feasible set of sub-channel patterns for user k can be summarized in the following matrix. 

01000100101 
00100110111 
00010011111 

00001001011 
where each row corresponds to the sub-channel index, and each column corresponds to the feasible sub- 
channel allocation pattern. Note that all the K users have the same allocation patterns matrix. We define 
a KJ indicator vector i = [ii, \k\ T where i& = [ik,i, ■-, ik,j] T , and where J is the total number of 
allocation patterns. Each entry i^j € {0, 1} which indicates whether a sub-channel pattern j is allocated 
to a user k or not. Since a single sub-channel pattern can be allocated to each user, maximizing the users' 
sum-utility is equivalent to maximizing the sum-utility of all users over all sub-channel allocation patterns 
such that each user is assigned a single pattern while respecting the exclusive sub-channel allocation 
constraint. Based on this analysis we have the following lemma. 

Lemma 3.1: The sum-utility maximization problem can be written as the following binary-integer 
programming problem: 



max < 



k=l j=l 



(6) 



K J 

s.t. ^2^2ik,jA k nJ = 1, Vn &) 
k=yj=i 

i fcii 6{0,l}, Vfc,j ©) 
where Uk,j{^ e Jj), a monotonically increasing function of the effective SNR 7^ is the utility of user k 

when allocation pattern j is chosen, and A\ j denotes the element of matrix A k corresponding to nth row 

and jth column. 

Proof: The proof is simple and follows from the following illustration. The effective SNR 7^/ of 



user k for pattern j is defined as: 

eff 



\- 1 



peak fjT 



1 



(7) 



where A^j is the number of sub-channels allocated to user k when allocation pattern j is chosen. The 
constraint ensures the exclusive sub-channel allocation i.e., any two sub-channel patterns allocated to 
two different users must not have any sub-channel in common. The constraint d6j)) means that at most one 
allocation pattern is chosen for each user. The per-user total power, the per sub-channel peak power and 
the allocated sub-channels power equality constraints are all implicitly accommodated in J^ff- ■ 

B. Joint Adaptive Modulation and Sum-Cost Minimization (JAMSCmin) 

1) JAMSCmin Problem Formulation: We now formulate the joint resource allocation and adaptive 
modulation problem. The objective is to allocate powers and sub-channels, and to choose the modulation 
scheme for each user in order to minimize the sum-cost while satisfying the the target data rate constraint 
of all the users (i.e., R^,\/k). For a modulation m G M to be chosen, the effective SNR of the user should 
not be less than a minimum value T^ that guarantees a target Block Error Rate (BLER) at the receiver. 
In addition, the power on all the sub-channels allocated to a user should be equal [2]. In the uplink the 
users terminals are more sensitive to transmit power due to their batteries's power limitations. Therefore, 
we introduce in the JAMSCmin formulation a user's cost which is function of its transmit power and has 
to be minimized. We define the following cost function for each user k 

C k (Pr x ,Pk) = - exp [Pr x - Pk] (8) 
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where p™ ax is the maximum power a user can transmit, and P k = YlneN * s tne sum °^ P owers 
transmitted by user k on its allocated set of sub-channels A4- The cost function is monotonically increasing 
in Pk whereas it is monotonically decreasing in p™ ax . With this choice of cost function, the JAMSCmin 
problem will not only minimize the sum-power of the users but will also ensure that each user's transmit 
power is minimized in accordance to its p™ ax level. In other words, a user with small p™ ax will transmit 
small power compared to another user with high p™ ax , and vice versa. The joint optimization problem 
can now be formulated as follows 

K 

min ^C k (Pr x ,Pk) (9) 

k=l 

s.t. R k >R%,Vk 

Pk,n = Pk,i,Vk,n,l 
Ik > T* m ,Vk,m 
\M k nM\ = l,Vfc 
A/j; n A/} = 0, Vfc ^ j 
nn( [j Mj) = | n £ {ni,ni + 1, ...,n 2 - l,n 2 } Lvfc 

where R k is the /cth user achieved data rate, M. k is a non-empty one element set that contains the modulation 
chosen for kth user; and where A4, and n 2 are the same as defined for SUmax problem. The fourth 
constraint reflects that a single modulation scheme is chosen for each user from the set M. In addition to 
its inherent difficulty due its combinatorial nature as explained for the SUmax problem, the joint adaptive 
modulation in addition to resource allocation renders the optimization problem © far more difficult to be 
solved. 

We now formulate this joint optimization problem as an equivalent BIP problem in the following. 

2) Equivalent BIP for JAMSCmin Problem: The sub-channel allocation patterns matrix is exactly the 
same as that for the SUmax problem. However, as the JAMSCmin problem considers joint adaptive 
modulation and resource allocation, we integrate the modulation selection into the sub-channel allocation 
patterns matrix. Since the number of sub-channels needed for transmitting a certain number of bits depends 
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B 



on the modulation scheme used, we refine the feasible allocation pattern matrix according to the modulation 
schemes. For example, the minimum number of sub-channels/TTI needed for = 140kbps is 3, 2 and 
1 for QPSK, 16QAM and 64QAM respectively. We recall that a TTI = 0.5msec, and each sub-channel 
contains 12 sub-carriers. Thus, for example, for total number of sub-channels N = 4, the kth user's feasible 
matrix of sub-channels allocation patterns for QPSK can be written as 

11111111101 
11111111111 
11111111111 
11111111011 
where the subscript m in corresponds to the modulation index. This matrix reflects that for the given 
RE, the number of sub-channels allocated to user k should not be less than 3 if QPSK is chosen. The 
same approach can be used to define kth user's sub-channels allocation patterns matrices for 16QAM 
and 64QAM. Depending upon their target data rates, the sub-channels allocation patterns matrices can be 
defined for all users on all modulation schemes. We define a KM J indicator vector £ = [£\ t i, ■■■,£k,m] T 
where £j cm = [4,mi> £km,j] T > an d where J is the total number of columns in the allocation pattern 
matrices. Each entry l^ m j € {0, 1} which indicates whether a sub-channel pattern j corresponding to 
pattern allocation matrix is chosen or not. Since a single sub-channel pattern and a single modulation 
scheme can be chosen for each user, minimizing the users' sum-cost is equivalent to minimizing the sum- 
cost of all users over all sub-channel allocation pattern matrices such that each user is assigned a single 
pattern and a single modulation scheme while respecting the exclusive sub-channel allocation constraint. 
Lemma 3.2: The joint resource allocation and adaptive modulation problem can be written as the 



following BIP problem: 1 K M J 

min \ g(£) = Y Y Y ^mjCkj^Pr^Pk,^) 

k=l m=l j=l 



KM.) 

k=l m=l j=l 
M J 

y ] y @k,m,j = ljVA; 

m=l j=l 

4,mJ G {0, l},Vfe,m, j 



1, Vn 



(10) 

CSt) 

Gob) 
Cot) 
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where n • denotes the element of matrix corresponding to nth row and jth column, P kj 



,m,3 — 



fi^k mj^k^m) * s me P ower transmitted by user k when jth sub-channels allocation pattern corre- 
sponding to B^ is chosen, and C k<m j(P% iax , Pk, m ,j) = ~ exp [P^ ax - Pk,m,j}- 



Proof: The transmit power Pk,m,j is a function of R^, T^, and the effective SNR J^ ^j of user k for 
jth pattern of B^. Let Pk,m,n be the power for user k on sub-channel n when modulation m is chosen, 
then jl'^j is g iven by x 




N k , mJ 2^n&N k , m , l+P k , m ,„G k , 



where Nk, m ,j with cardinality iVjt m is the set of sub-channels allocated to user k when jth pattern from 




is chosen. The power allocation values Pk, m ,j' s are obtained prior to resource allocation by solving 



the following equations: 




) 




0,Vk,m,j 



(12) 



which are obtained by setting j k 



eff _ 
k,m,j 



T* m and P k . 




,3 



^ and hence the per user minimum SNR and 



the allocated sub-channels powers equality constraint are implicitly accommodated in Pk, m ,j- The per-user 
target data rate constraint is already implicitly accommodated in the definition of allocation patterns and 
hence in the calculation of P k ,m,j- The constraint (TTOb ) reflects the mutual exclusivity restriction on the 
sub-channels allocation and constraint (TTOb) means that at most one allocation pattern and one modulation 



We recall that the formulation of the problems as equivalent binary-integer programs is an intermediate 
step towards their solution. Although the BIP problems may look simple compared to the primal problem 
but unfortunately, their solutions are exponentially complex due to their combinatorial nature. A similar 
binary-integer programming solution was proposed for weighted-sum rate maximization problem in lTT6l 
but as mentioned before it is exponentially complex which is not practical. In the following section, we 
propose a polynomial-complexity framework for the solution of both the above problems that is inspired 
from the canonical dual transformation method. The main idea of our proposed approach is to transform 
each binary-integer programming problem into a canonical dual problem in the continuous space whose 
solution is identical to the corresponding binary integer program under certain conditions. 



scheme is chosen for each user. 
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IV. Canonical Dual Approach for Solving the BIP Problems 



Under certain constraints/conditions, the canonical duality theory [18] can be used to reformulate some 
non-convex/non-smooth constrained problem into certain convex/smooth canonical dual problems with 
perfect primal/dual relationship. However, this theory does not provide any general strategy for the solution 
of non-convex/non-smooth problems. The constraints under which the canonical dual problem could be 
perfectly dual to its primal problem is purely dependent on the nature of the primal problem under 
consideration and should be studied for each specific problem anew. This theory comprises of canonical 
dual transformation, an associated complementary-dual principle, and an associated duality theory. The 
canonical dual transformation can be used to convert the non-smooth problem into a smooth canonical 
dual problem; the complementary-dual principle can be used to study the relationship between the primal 
and its canonical dual problems; and the associated duality theory can help to identify both local and global 
extrema. Comprehensive details about this theory, and its application to an unconstrained 0-1 quadratic 
programming problems can be found in ||l"8l , and ||T9l respectively. Due to the presence of additional 
constraints, our problems are far more difficult compared to that described in [19]. 

By using the aforementioned theory, we transform each of the SUmax and JAMSCmin primal problems 
into a continuous space canonical dual problem in the following. We then study the optimality conditions, 
and prove that under these conditions, the solution of each canonical dual problem is identical to that of 
the corresponding primal problem. 

A. Canonical Dual Problem and Optimality Conditions for SUmax Problem 

The objective function, V(\) in problem © is a real valued linear function defined on I„ = i C M. KxJ 
with feasible space defined by 

K J J \ 

l f = liel a cR KxJ | J2J2 i kd A n,j = 1 ^ E^i = 1 ' VA; ; **je{0,l}Vfc,4 (13) 

k=lj=l j=l ) 

We start our development by introducing new constraints ikj{ik,j — 1) = 0,V/c,j which means that 
any ik,j can only take an integer value from the set {0, 1}. This approach is used for the solution of a 
0-1 quadratic programming problem in [19|. However, the problem considered in [19] is a simple uncon- 
strained 0-1 quadratic programming problem while our problem is combinatorial in nature with additional 
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constraints. In other words, in addition to the binary-integer constraint on ijt j's, we have the mutual exclu- 
sivity restriction on the allocation of sub-channel patterns (i.e., {i^j x iij = 0\k ^ l,\/k, I £ {1, K}}), 
and the mutual exclusivity constraint on the sub-channel allocation i.e., X^feLi X^f=i ^k,j^n j = l,Vn. 
Furthermore, at most one sub-channel pattern can be allocated to a user i.e., ^/=i**,j = IjVA;. Note 
that the mutual exclusivity restriction on the sub-channel patterns allocation is accommodated implicitly 
in the formulation of the primal problem and does not show up explicitly. We temporarily relax the new 
constraints ik,j(ik,j ~ 1) = 0,Vfc,j, and the equality constraints (|6H6])) to inequalities and transform the 
primal problem with these inequality constraints into continuous domain canonical dual problem. We will 
then solve the canonical dual problem in the continuous space and chose the solution which lies in If as 
defined by (fT3l ). Furthermore, for our convenience, we reformulate our primal problem as an equivalent 
minimization problem. The primal problem with these inequality constraints can now be written as follows. 



mm < 
i 



K J 



ff^-EE* (14) 



k=l j=l 
v K J 



3=1 

ik,j {ik,j - 1) < 0, Vfc, j 



where Ukj is used to denote Ukj("jt /) and will be used in the remainder of the paper. 



i kJ £ {0, 1}, Vk,j 

■j is useu lu uenuie UkjVlk j ^ 
The temporary relaxation of the constraints to inequalities is needed for developing the canonical dual 
framework. We prove later that the solution of the canonical dual problem achieves the binary-integer 
constraints i.e., ikj (ik,j — 1) = 0,Vk,j and all the other constraints with equality. As a first step towards 
its transformation into a canonical dual problem, we relax the primal problem ( lTT~8l . [19]). To this end, 
we define the so-called canonical geometrical operator x = A(i) for the above primal problem as follows: 

x = A(i) = (e, A, p) : R KJ -> R N x R K x R KJ (15) 

which is a vector- valued mapping and where p = [if(ii — 1), — 1)] T is a KJ-vector with ij[ (ij. — 

1) = [*fc,l(*jfc,i - 1), -,ik,j(ik,J ~ 1)] T > ^ = fe/=i*i,j ~ x ) ' (E/=i^J - is a K-vector and 
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[Y,k=l E/=i ikj^ij - l) , (EfcLi E/=i '/.-., -^v.j - l) is an N- vector. Let Xa be a convex 
subset ofx = M iV xf p x R KJ defined as follows 

Xa = {x = (e,\,p) eR N xR K xR KJ | e<0,A<0,p<0} (16) 

We introduce an indicator function V : x ~ * ^ U {+°o}, denned as 

if x G Xa, 



y(x) = < 



(17) 

+oo otherwise. 



Thus, the inequality constraints in the primal problem ([Pil l can now be relaxed by the indicator function 
V(x), and the primal problem can be written in the following canonical form [19]: 

min\v(A(i))-J2J2 i kd U k,j I ikj€{0,l}Vk,j\ (18) 
{ fc=ij=i J 

We now define the canonical dual variables and the canonical conjugate function associated to the indicator 
function in order to proceed with the transformation of the primal problem into canonical dual. Since V(x) 
is convex, lower semi-continuous on x> the canonical dual variable x* G x* = X = ~^> N x ^ A x R KJ is 
defined as: 



(e*,\*,p*) if e* > € R N ,X* > G R K , p* > G 
x*£dV(x) = { ~ (19) 

otherwise. 

By the Legendre-Fenchel transformation, the canonical super-conjugate function of V(x) is defined by 

V $ (x*) = sup{x T x*-y(x)} =supsupsup{e T e* + A T A* + p T p*} 
xex c<OA<Op<o 

if e* > 0, A* > 0, p* > 0, 

(20) 

+oo otherwise. 
The effective domain of V"(x) is given by 

x * a = {(e*,\*,p*) (£R N xR K xR KJ | e* > G R N , A* > G R K , p* > G R KJ } (21) 

Since both V(x) and 0(x) are convex, lower semi-continuous, the Fenchel sup-duality relations 

x* G dV(x) x G dV*(x*) V(x) + y»(x*) = x T x* (22) 

hold on x x X* ■ The pair (x,x*) is called the extended / Legendre canonical dual pair on \ x X*, an d 
the functions V(x) and V"(x) are called canonical functions lfT8l . The optimal solution of our primal 
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problem can be obtained if and only if x = If G Xa, i.e., along with the satisfaction of the binary-integer 
constraints, all the other constraints must be achieved with equality. Thus, we need to study the conditions 
under which the canonical dual variables x* G x* a can ensure that x = If £ \a- By the definition of 
sub-differential, the canonical sup-duality relations (l22l are equivalent to the following: 

x < 0, x* > 0, x T x* = (23) 

From the complementarity condition x T x* = 0, for x* > 0, we have x = (i.e., e = 0, A = 0, p = 0) and 
consequently x = If G Xa- This means that for x* > 0, all the constraints of the primal problem (TT4T > are 
achieved by equality (with i^j G {0, 1}, V/c, j which comes from p = 0). Thus, the dual feasible space for 
the primal problem is an open positive cone defined by 

X ; = {(e*,X*,p*)€x*a I e* >0,A* >0,p* >0} (24) 

The so-called total complementarity function (see lTT8l , |[T9l for definition), E(i, x*) : x x X\ ~ * K associated 
with the primal problem (fT4l l can be defined as follows. 

K J 

E(i,x*) = A(i) T x* - V\x*) -YtYt^jU** (25) 

fc=l j=l 

which is obtained by replacing V(A(i)) in (fT8l) by A(i) T x* — V"(x*) from Fenchel sup-duality relations 



d22l) . From the definition of A(i) and V"(x*), the total complementarity function takes the form: 

K J r / N \ \ N K 

S(i, £*, X*,p*) =J2J2\ P*kA,j +[ X *k- Pt,j ~ U k,j + € n A n,J ik >3 \ ~ Yl ^ ~ Yl A fc ^ 
k=lj=l I V n=l / J n=l fc=l 

Similar to [19], the canonical dual function f d (e* , \* , p*) associated to our primal problem for a given 
(e*, A*, p*) G xj can be defined as 

f d (e*,A*,p*) = sta{H(i,e*,A*,p*) | i G I a } (27) 

where sta{/(x)} stands for finding the stationary points of f(x). The complementarity function is a 
quadratic function of i G I a , and has therefore a unique stationary point with respect to it for a given 
(e*, A*, p*) G x*. The stationary points of 3(i, e*, A*, p*) over i G I a occurs at i(x*) with 



*fcj( x *) 



^ \U k ,3 + pi, -K~Y e n A n,l) » yk >i ™ 



16 

Replacing ikj by i^jix*) in (l26l ). we have 

f (e , a , P*) = - i EE — ^ — - E < - E A * w 

k=l j=l I ' k j I n=l fc=l 

which is a concave function in xj- The canonical dual problem associated with the primal problem ([141 ) 
can now be formulated as follows 

ext{fV,A*,p*) I (e*,A*,p*) exj} (30) 

where the notation ext{/(x)} stands for finding the extremum values of f(x). 

We have the following canonical duality theorem (Complementary-Dual Principle) on the perfect dual 
relationship between the primal and its corresponding canonical dual problem. 

Theorem 4.1: If (e*, A*,p*) G xj is the stationary point of f d (e*, A*, p*), such that 

i= [i 1A ,...,i K ,jf with 1^3 = —^- lUkj + plj-Jl-^relAnj) ,Vfc,j (31) 

\ n=l / 

is the KKT point of the primal problem, and 

f(i) = f d (e*,r,p*). (32) 

then the canonical dual problem (l30l is perfectly dual to the primal problem ©. 

Proof: The proof can be obtained directly from the proof of Theorem 1 given in 11191 . ■ 
The above theorem shows that the binary-integer programming problem (fTOl is converted into a dual 
problem in continuous domain which is perfectly dual to it. Furthermore, the KKT point of the dual 
problem provides the KKT point for the primal problem. However, as the KKT conditions are necessary 
but not sufficient for optimality in general, we need some additional information on the global optimality. 
Based on the properties of the primal and dual problems, we have the following theorem on the global 
optimality conditions. 

Theorem 4.2: If (e*,A ,p*) G x|, then i defined by (HQ is a global minimizer of f(i) over Xf and 
(e*, A*,p*) is a global maximizer of f d (e*, A*,p*) over x|, and 

f(i) = minf(i) = max f d (e*, A*, p*) = f d (e*, A*, p*). (33) 

(e*,A*,p*)GxJ 

Proof: See Appendix A. ■ 
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B. Canonical Dual Problem and Optimality Conditions for JAMSCmin Problem 

The difference between the SUmax and the JAMSCmin problems lies only in their formulation. The 
procedure and steps of canonical dual transformation for JAMSCmin problem are the same as that followed 
for the SUmax problem, and therefore, we will not repeat them in this paper. The canonical dual problem 
associated to the JAMSCmin primal problem (fTOt can be obtained as given by 

max{ g d (r,/Af?*) i (e^*,Q*)^y;} m 

where yZ is the associated dual feasible space defined as 

y* = {(£*, q*) £R N x R K x R KMJ | £* > G R N , fx* > G R K , Q* > € R KMJ } (35) 

and g d (£*, n* , Q*) : 1^ x R K x R KMJ — >■ R is the corresponding canonical dual function defined as 
follows: ( / „ A r , \ 2 ^ 

g d (^ ,m ,q ) = -^2^2^2^\- ~* — \-Zs&-Zsi J * (36) 

k=l m=l j=l I y fc,m,j n= l fc =1 

which is a concave function on yt, and where Cfc m j is used to denote Ck tm j(P™ ax , Pk, m ,j)- Moreover, 
the results on the primal/dual relationship (perfect duality) and the global optimality conditions can be 
obtained by a similar procedure followed for SUmax problem (i.e., Theorem 4.1, and Theorem 4.2). 

Based on the above mathematical analysis, we provide resource allocation (with joint adaptive modulation 
for JAMSCmin) algorithms in the following section. An adaptive modulation scheme for SUmax problem is 
also proposed since unlike the JAMSCmin problem, it does not capture the adaptive modulation implicitly 
in the problem formulation. The proposed adaptive modulation is based on the powers and sub-channels 
allocated to each user by the proposed resource allocation algorithm. 

V. Resource Allocation and Adaptive Modulation Algorithms 

A. Resource Allocation Algorithm for SUmax 

The proposed algorithm is based on the solution of canonical dual problem which according to theorem 
4.2 provides the optimal solution to the primal problem if the given global optimality conditions are met. 
Since the dual problem is a concave maximization problem over \\ > it i s necessary and sufficient to solve 



18 



the following system of equations for finding the optimal solution iBTTl . 

B rd K J ( / N \ 1 



d\* k ^ 2p* k . 

K J = l K >J \ 71=1 



di d _ /^-A^-Eti^,; 



2 

-1 = 0, VkJ (39) 



We propose a sub-gradient based iterative algorithm for the above system of non-linear equations that 
is equivalent to solving f d (e*,A*,p*) using gradient-decent method [21 |. The interest of using the sub- 
gradient method is its ability to use the decomposition technique that allows to simplify the solution by 
using a distributed method. The iterative algorithm is given in Table U where each of q, s and t denotes the 
iteration number and j3 p * , j3\* and /? e » denote the step sizes for the sub-gradient update. For an appropriate 
step size, the sub-gradient method is always guaranteed to converge [21]. The algorithm starts by initializing 
the variables. Then, for the given e*^ and X*^°\ the solution to the set of equations d39l i.e., p*te) is 
obtained in step 1. The operation 



(40) 



' • sgn (/;;;. /" ; 'J // if (p*^ 1 ) +(3 p ,C {q - 1) ) = 0,Vk,j 
p* (9) = p*^- 1 ) + (3 p ,C {q - 1] otherwise. 



in step 1 is the projection of p* onto the space x* P ' = {p* € M. KJ \p* ^ 0}, since the canonical dual 
objective function is not defined at p* = 0. In d40b , sgn stands for sign/signum function and < ry << 1. 
According to the above projection, if the updated value of p* in the current iteration occurs to be zero, 
it is projected to the negative domain if its value was positive in the previous iteration, and vice versa. 
This projection has no impact on the convergence, since the sign of p* does not change the direction of 
the gradient (see equation (|39l). Step 2 finds A*^** that solves equations' set (1381 ) for the given e*(°) and 
p* ^ . These values of p* ^ and A* ^ are then used to solve the set of equations (f3VT > by updating e* ^ to 
e* w in step 3. Step 4 checks whether ||^| < 5 for p* (g) , A* (s) and the updated e* (t) where 5 -> is the 
stopping criterion for sub-gradient update. If | J-j— | > 5, steps 2 through 4 are repeated until both [ j^- \ < 5 
and |gr| < 5. In step 6, C (g) is recomputed for p*^\ and the updated A* (s) and e*W. If |C (g) | < 5, the 
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algorithm is stopped otherwise steps 1 through 6 are repeated until convergence. The resource allocation 
vector i is then obtained from the dual optimal solution (e*, A*,p*) in step 8. 

TABLE I 

Resource Allocation Algorithm 

Initialize (e* (0) , A* (0) , p* (0) ) G X J 

1. Compute C (9) = f^lp-w- If IC (9) I < S, go to step 2. 

.set p*« +1 > ^n Xp , (V w +/W 9) )- 

• Set q «— q + 1, and repeat step 1. 

2. Compute r) {s) = |^| A , (s) . If |f? (s) < 8, go to step 3. 

.Set A*< a+1 >^ (a*W+/W s) )- 

• Set s <— s + 1, and repeat step 2. 

3. Compute u (t) = f£| e „(t). If |w (t) | < 5, go to step 4. 

.Sete*( t+1 »<-(e*W+/3 £ . W (*)). 

• Set t t + 1, and repeat step 3. 

4. Recompute T7 (s) = J^r A ,( s >- 

5. Repeat steps 2 through 4 until |?7 (s) | < S, and |w (t) | < 6 

6. Recompute C W = f£| p .w 

7. Repeat steps 1 through 6 until |C <9) | < 5, |ry (s, | < 5, and |t> (t, | < S 

8. Compute i according to (31\ . 



1) Adaptive Modulation Scheme for SUmax: By knowing perfectly the effective SNR of each user from 
the powers and sub-channels allocation performed according to the previous subsection, we propose an 
adaptive modulation scheme in this subsection. The proposed adaptive modulation scheme is based on 
the criterion of target Target Block Error Rate (BLER) at the receiver used for the JAMSCmin problem. 
According to this approach, for a modulation m G M to be chosen, the effective SNR of the user should 
not be less than a minimum value Y* m that guarantees a target BLER at the receiver. Since the effective 
SNR of users are perfectly known from the the powers and sub-channels allocation performed according 
to the previous subsection, we adopt the modulation for each user which maximizes its individual utility. 
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Thus, depending upon 7^, the efficient modulation for user k is determined as follows: 

m*(k) = arg min {(7^ - T* m )\ r , e„| (41) 

Note that the above approach is similar in spirit to the approach used in l20l where adaptive modulation 
in OFDM system is considered and an efficient constellation is chosen for each sub-channel. 

B. Joint Adaptive Modulation and Resource Allocation Algorithm for JAMSCmin 

The dual function g d (£* , p*, g*) is a concave function over (£*, p*, g*) £ yi. Thus, the corresponding 
dual problem is a concave maximization problem over y% where the joint adaptive modulation and resource 
allocation can be obtained by solving the following set of equations: 

3 d K M J f / N \ ^ 

= £ £ £ t^— ei ind - c k , m>j - ni - J2 CB k m , n , 3 Bt^ -1 = 0, 

k=lm=lj=l { Ze k,m,j \ n=l / J 



o d M J ( , / N 

Og 



£ £ 1 On* ( ~~ Ck,m,j ~ lA ~ £ Cn^m,n,j J \ ~ 1 — 0; ^ 



m=l j = l K. K >3 
9g d ( ~Ck,m,j ~ P% ~ En=l e n B m,n,j \ 



(43) 



n=l 



1 = 0, Vk,m,j (44) 



^@k,m,j V @k,m,j J 

A similar procedure of sub-gradient is proposed where an iterative algorithm can be derived that is 
similar in spirit to that derived for the SUmax problem. Since it uses a similar procedure and has a similar 
sequence of steps as that for the algorithm given in Table I, the latter can be adopted to the JAMSCmin 
problem, and we do not reproduce it in this paper. 

C. Complexity of the algorithm 

1) Complexity of the algorithm for SUmax problem: In each iteration for p* , we compute KJ variables. 
The number of variables computed in each iteration for A* is K and that for e* is N. Assume that the 
number of iterations required for optimal p* , A* and e* are l p -, 1\* and I e * respectively, then the algorithm 

has an overall complexity of 0(I p *KJ + l\*K + I e *N). 

2) Complexity of the algorithm for JAMSCmin problem: The complexity of the proposed algorithm 

adopted to the JAMSCmin problem is <D(I e *KMJ + I^K + I^N) where I e *, 1^ and 1^ are the 
numbers of iterations needed for finding the optimal values of KM J variables g* , K variables p* and 
the N variables £* respectively. 



21 

D. On the Optimality of the Algorithm 

The canonical dual problem is a concave maximization problem over xt> the proposed algorithm is then 
surely optimal if (e*,A*,p*) € xj- However, if (e*,A*,p*) is not inside the positive cone xj, then the 
canonical problem is not guaranteed to be concave. Consequently, the proposed algorithm may not find the 
optimal solution. From our simulation results, we have observed that for moderate number of sub-channels 
the proposed algorithm works well, and the canonical dual solution is very close to the optimal solution. 

In this subsection, we analyze the gap between the optimal solution and the solution obtained by using 
our proposed sub-gradient based algorithm. We perform the analysis for SUmax problem which is equally 
applicable to the JAMSCmin problem, and we will not repeat it in this paper. We start the analysis by 
introducing a modified problem whose optimal solution is not necessary and will not replace our actual 
problem but is used only to study the optimality gap of our proposed algorithm. In our analysis, first 
we find the solution of the modified problem (which is a stationary point and may not be necessarily the 
optimal solution of this modified problem). Then, we show in Theorem [5J] that there exist a primal problem 
with a slightly different values of the utilities Ukj's whose optimal solution is equal to the solution of 
this modified problem. Finally, in Corollary 15. II we show that under certain conditions, the solution of the 
canonical dual problem obtained using the algorithm in Table U provides solution to the primal problem 

which is very close to optimal solution. Let us consider the following modified problem 

(VI) : max i d (e\\\p*) (45) 

s.t. e* > c g5t) 
A* > d (@3>) 

where (c, d) € (M+ , ). We solve this problem using the standard Lagrangian technique. Let (e*, A*, p*) 
be the obtained solution. The corresponding Lagrangian can be defined as 

L(e*, X*,p*,cr e " , cr A * , aP') = f d (e*,\*,p*) - (e* T - c T )a e ' - cr A * (A* T - d T ) (46) 

where (er e *,cr A *) G (M^K^) are the Lagrange multipliers associated to the constraints (I45k-l45b) respec- 
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tively. The corresponding KKT conditions are: 

K J ( -, / N 



dL 
oe* 



dL 



E E [ 2^- + ^ - A * - E «i J <i | = 1 + > Vn ^ 47 ) 



j=l ^ K >J \ n=l 



9L =► f ^' ^ ) _ 1 = , VA-.j (49, 



The above equation can be solved using the sub-gradient based algorithm in Table U Moreover, in order 
to ensure that the solution of (|49l ) is obtained for positive p*, we can use the following projection in the 
update of p*: 

= argmin p . 3 . 6x j U^T 1 ) - p y\ if <& fcj .(i-D < ,Vfc, ' 



p * (q) _ ® kj (i- 1 ) otherwise. 



(50) 



where = p*(<?— !) + /3 p » M±.\ p , {q _ 1} denotes the sub-gradient update, and where /3 P » is the step size. 

The above projection ensures the positivity of p* . 

Theorem 5.1: For U^j = U^j — 2&k,jPkj w i tri @k,j = { — 1)0) !}• ^iJ! there exists a primal problem 
f(i) with utilities Ukj replaced for Ukj that can be solved optimally using the algorithm in Table I. The 
solution (e*,A*,p*) of f(i) obtained using Table I is equal to the solution of the modified problem (l45l l. 

Proof: See Appendix B for the proof. ■ 
Moreover, we have the following result which is the corollary of Theorem 5.1. 

Corollary 5.1: If p* k ■ « fT^j, Vk,j; then the solution of the canonical dual problem obtained using 
the sub-gradient based algorithm (Table I) provides a solution to the primal problem which is very close 
to the optimal solution. 

Proof: See Appendix C for the proof. ■ 
1 ) Analysis of the algorithm 's results for N — > oo: It can be seen from the KKT equation ( f38l) that 

U k,j ~ Afe - En=i £*n A n,j = P*k,j when a P attern 3 is allocated to user k, and U Kj - A* - Y,n=i e* n A k n j = 
—p* k ■, otherwise. When the number of sub-channel is very high, there are several patterns that have nearly 
equal utilities £/& /s for each user. This is due to the fact that for high number of sub-channels, the 
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per sub-channel utility will be very small, and since the difference of sub-channels in the patterns with 
high number of sub-channels will be less, their utilities will have very small difference. Furthermore, the 
difference between the summation term J2n=i e n A n j f° r sever al patterns of user k will be very small. 
This means that the term Ukj — A|, — J2n=i e n A n j f° r several patterns of user k will be nearly equal, as 
A£ is the same for all the patterns of that user. 

Let us assume that a pattern j is allocated to user k. Consequently, Ukj — A|, — J2n=i e n A n j = P*k j 
whereas Uk,j' — X* k — Yln=i e n A tj' = ~~ Pkj' ^ or a ^ ^ 3- Moreover, in view of the above discussion, 
the difference between U k j - X* k - Yln=i e n A n,j and u k,j' ~ K~ En=i e n A n,j> wiU be ver y sma11 
for j and j' with high number of sub-channels. Thus, being equal to Ukj — XI — J2n=i e n A nj and 
—Uk,j> + A£ + Yln=i e n A nj' respectively, p* k ■ and p\ ■, will both be very small compared to Ukj and 
Ukji respectively. 

VI. Numerical Analysis and Discussion 

We consider a system with 5MHz of bandwidth (i.e. LTE) divided into TV = 25 sub-channels each 
having a bandwidth of 180kHz. We assume that K = 10 uniformly distributed users are simultaneously 
active in the cell. The scenario assumed is urban canyon macro which exists in dense urban areas served 
by macro-cells. A frequency selective Rayleigh fading channel is simulated where the channel gain has a 
small-scale Rayleigh fading component and a large-scale path loss and shadowing component. Path losses 
are calculated according to Cost-Hata Model and shadow fading is log-normally distributed with a standard 
deviation of 8dBs. The power spectral density of noise is assumed to be -174dBm/Hz. 

A. Sum-utility maximization 

In simulations, we assume that the utility of the user is equal to its weighted rate where the rate is 
defined by Shannon's formula. In other words, the SUmax problem is equivalent to weighted-sum rate 
maximization. Fig. 1 plots the empirical cumulative distribution function (CDF) of sum-utility for different 
resource allocation algorithms. The figure illustrates the comparison of the CDF's corresponding to our 
proposed algorithm, both the binary-integer programming solution and the greedy algorithm proposed in 
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lfl6l . and the round robin scheme in which an equal number of consecutive sub-channels are allocated 
to each user in turn. The figure shows that although the greedy algorithm proposed by Wong et al. is 
efficient in comparison to the round robin scheme, its performance is far away from the proposed solution. 
Moreover, it can be seen from the figure that the results of the proposed algorithm are very close to that 
obtained by solving the binary-integer program which is the optimal solution. 

B. Joint Adaptive Modulation and Resource Allocation 

The minimum effective SNR for each modulation T* n that ensures a target Block Error Rate at the 
receiver is determined from the link-level performance curves (e.g., see ll22lD . Fig. 2 displays the empirical 
CDF of sum-cost for different algorithms when sum-cost minimization based resource allocation (RA) is 
performed joint with and without adaptive modulation (AM). The figure illustrates the comparison of the 
CDF's corresponding to our proposed resource allocation algorithm when joint AM and RA is performed 
and when RA is performed without AM, the binary-integer programming (BIP) based solution adopted to 
joint AM and RA problem, and RA without AM, and the round robin scheme in which an equal number 
of consecutive sub-channels are allocated to each user in turn and minimum possible power is allocated to 
the users while ensuring their target data rates. The round robin scheme is used as a baseline scheme for 
comparison. The RA without AM scenario considers 16QAM as the modulation scheme. The proposed 
RA with fixed modulation outperforms the round robin scheme which is not unexpected. The figure shows 
that the joint AM and RA results in a significant performance improvement over the RA without AM. 
The performance of the proposed algorithm can be depicted from the fact that the results of the proposed 
algorithm nearly overlap with that of the BIP based solution both for joint AM and RA, and RA with 
fixed modulation scheme. We recall that the BIP based solution is the optimal solution. 

VII. Conclusion 

This paper studies resource allocation and adaptive modulation in uplink SC-FDMA systems. Sum-utility 
maximization, and joint adaptive modulation and sum-cost minimization problems are considered whose 
optimal solutions are exponentially complex in general. A polynomial-complexity optimization framework 
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that is inspired from the recently developed canonical duality theory is derived for the solution of both the 
problems. Based on the resource allocation performed by the proposed framework, an adaptive modulation 
scheme is also proposed for the sum-utility maximization problem that determines the best constellation 
for each user. The optimization problems are first formulated as binary-integer programming problems 
and then, each binary-integer problem is transformed into a canonical dual problem in the continuous 
space which is a concave maximization problem. The transformation of the problem in continuous space 
significantly improves the performance of the system in terms of complexity. The proposed continuous space 
optimization framework has a polynomial complexity that is a significant improvement over exponential 
complexity. It is proved analytically that that under certain conditions, the solution of the canonical dual 
problem is identical to the solution of the primal problem. However, if the dual solution does not satisfy 
these conditions then the optimality can not be guaranteed. Therefore, some bounds on the sub-optimality 
of the proposed framework when these conditions are not satisfied are also explored. The performance 
of the proposed canonical dual framework is assessed by comparing it with the existing algorithms in 
the literature. The numerical results show that the proposed framework provides integer solution to each 
problem which is very close to optimal. 



A. Proof of Theorem 4.2 

The total complementarity function S(i, e*, A*, p*) is convex in i and concave (linear) in e*, A* and p*. 
Therefore, the stationary point (i, e*, A*, p*) is a saddle point of S(i, e*, A*, p*). Furthermore, i d (e*, A*, p*) 
is defined by S(i, e*, A*, p*) with i being a stationary point of S(i, e*, A*, p*) with respect to i G l a . 
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Consequently, f d (e*,A*,p* 



) is concave on \\ an d the KKT point (e*,A*,p*) G \\ must be its global 



maximizer. Thus, by the saddle mini-max theorem: 



i d (e*. A*, p*) = maxmaxmaxf d (e*. A* 



) 



e*>0 A*>0p*>0 



= max max max min e* 

e*>0 A*>0p*>0 iex„ 




max max max mm 
e*>0 A*>0p*>0 ieJ a 



{f(i) + e r 6* + A T A* + pV} 
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max max mm < 

e*>0 A*>0 ieZ 



f(i) + e e* + \ A* + max < 

p*>0 



K J 
k=l j=l 



max mm 




f(i) + e T e* + E A M E " 1 

k=i \j=i 





mm J f(i) + max ^ < E E " 1 

a [ [n=l \fc=lj=l 

J 

s.t. ik,j(ik,j ~ 1) = 0,VA;, j; = 1,VA; 

i=i 

f J K J 

minf(i) s.t. { ik,j(ik,j ~ 1) = Q^V^iiE^ = ^^EE^J^i = 1)Vn 



j=i fc=i i=i 



minf(i) (51) 



Note that the linear programming 



max < 

p*>0 



K J 

fe=i j=i 



has a finite solution in the open domain xj if and only if ik j (ifc j — 1) = 0, VA;, j. By a similar argument, the 
solution of max A - >0 [zti K (l)/=i hj ~ l) } and max e , >0 {^ =1 4 (£f =1 £/ =1 tfcj^ " l) } 
leads to the last equation ([STJ. This shows that the KKT point (e*, A*, p*) maximizes f d (e*, A*, p*) over 
if and only if i is the global minimizer of f(i) over If. This completes the proof. 

B. Proof of Theorem 5.1 

Using sub-gradient method with projection defined by (l50l l ensures the positive solution of KKT equation 
d49l ) which implies that the corresponding ik j is binary integer. However, respecting the positivity constraint 
on A*, equation (l48l) can not ensure that a single sub-channel pattern is allocated to each user but 1 + a* 
number of patterns will be allocated to each user k. Similarly, ensuring that e* > 0, equation (l47l) means 
that a sub-channel can be allocated to more than one users. 

In the following, we discuss that we can find another approximate problem for which the above KKT 
equations not only provide binary integer solution but also ensure that a user will be assigned with a 
single sub-channel pattern and a sub-channel will be allocated to a single user. To this end, we proceed 
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as follows. The KKT equation (4% can be written as 

N 

Uk,j - e n A n,j = ± P*k,ji j (52) 

n=l 

We introduce KJ new variables flfcj's denned as follows 

{ {1, 0} if U kJ -XI- Zn =1 e* n A* nJ = -p%. 
Vk,j = < (53) 

[ {-1,0} iiUkj-Xl-ZLl^nK^+Pt, 

From the above definition of 9 k equations (l52l can be written as 

AT 

C/ fe)i -20 fc) X J -^-^e;4 ii = ±p^, Vfc.j (54) 

71=1 

Let Uk j = Uk j — 29k jPk then the above equations take the form: 

' 3 N 

Ukj -K-J2 e n A h = iPfcj, Vfc, j (55) 

n=l 

Although the utilities are changed from Ukj to fT^- = Uk j — 29kjP* k j> th e solution of the above equations 
provide integer solution to ifcj's. We now apply this change in utilities to the equations ( |47Tl48l ). The KKT 
equations (|48l ) can be written as 

Replacing [7^^- for JJfcj — 29k^p* k j, the above equations become: 

E { 2^- + P ^ " A * " E < A n)j } + E ^ = 1 + > vfe ^ 57 ) 

If there exist 0fcj's such that ^/=i ^fcj = ^fc*' then we have 

E { 2^7 + ^ " ^ " E < A n)j } = VA; ^ 58 ) 

This implies that there exist another problem with a different set of utilities for which the above solution 
ensures that a single pattern will be allocated to each user. By using a similar procedure for the KKT 
equations (|47T ). we get 

E E | 2^T + p h ~ ^ ~ E «i) <i | = 1, Vn (59) 

which enures that a sub-channel will be allocated to a single user at most when Y2k=l S/=i ^jAj j = <r^* , 
and the utilities are changed from Ukj to [T^j = U k> j — 29kjp k j. 

The above analysis shows that the solution of the problem VI, namely (e*, A , p*) that lies in the positive 
cone, is the solution of the above KKT equations <55I58I59|) . Moreover, the KTT equations (1551581591 ) give 
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the stationary point of a slightly modified problem P(e*, A*, p*) which is the canonical dual of a slightly 
modified primal problem with utilities Ukj = Ukj — 2dk,jP%j- Since the solution (e*, A*,p*) is positive, 
according to Theorems 4.1 and 4.2, the proposed sub-gradient based solution proposed in Table U optimally 
solves a corresponding primal problem with utilities Uk,/ S an d an objective function f(i). Note also that 
the canonical dual P(e*,A*,p*) is concave (since the KKT solution is in the positive cone). However, 
how far the solution of the modified problem will be from that of the primal problem © depends upon 
the values of p\ 's. 

C. Proof of Corollary 5.1 

If P*k,j « U kJ ,Vk,j, then U kJ ^ U kJ ,Vk,j, f(i) » f(i), and 

max f d « max f d (60) 

(e*,A*,p*) (e*,A*,p*) 

For p^, j << Uk j, Vk, j, the solution of the equations (I55I58I59I ) is very close to that of equations (1371381 
l39l ). Furthermore, the solution of (I55I58I591 I is the optimal solution of the corresponding primal problem 
with utilities Ukj (which is very close to the optimal solution of the primal problem with utilities Ukj). 
Consequently, the dual canonical problem obtained using the sub-gradient based algorithm (Table U) will 
provide solution to the primal problem which is very close to the optimal solution. This completes the 
proof. 
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Fig. 2. Empirical CDF of sum-cost 



